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Abstract 

Context. R Coronae Borealis stars (RCBs) are hydrogen-deficient and carbon-rich supergiant stars. They are very rare, as only ~ 50 
are actually known in our Galaxy. Interestingly, RCBs are strongly suspected to be the evolved merger product of two white dwarfs 
and could therefore be an important tool to understand Supemovae type la in the double degenerate scenario. Constraints on the 
spatial distribution and the formation rate of such stars are needed to picture their origin and test it in the context of actual population 
synthesis results. 

Aims. It is crucial to increase significantly the number of known RCBs. With an absolute magnitude My ~ -5 and a bright/hot cir- 
cumstellar shell made of amorphous carbon grains, RCBs are so distinctive that we should be able nowadays to find them everywhere 
in our Galaxy using publicly available catalogues. In the optical, the search is difficult as RCBs are known to undergo unpredictable 
photometric declines. However mono-epoch mid-infrared data can help us to discriminate RCBs among other dust-producing stars. 
The aim is to produce from the mid-infrared WISE and near-infrared 2MASS catalogues a new catalogue of reasonable size, enriched 
with RGB stars. 

Methods. Colour-Colour cuts used on all stars detected are the main selection criteria. The selection efficiency was monitored using 
the 52 known RCBs located in the sky area covered by the WISE first preliminary data release. 

Results. It has been found that selection cuts in mid-infrared colour-colour diagrams are a very efficient method of discriminating 
RCBs from other stars. An RGB enriched catalogue made of only 1602 stars, with a high detection efficiency of about 77%, was 
produced. Spectral energy distributions of 49 known RCBs and 5 known HdCs are also presented with estimates of their photosphere 
and circumstellar shell temperatures. 

Conclusions. The newly released WISE all sky catalogue has proven to be a valuable resource in finding RGB stars. Actual scenarios 
predict that between 100 and 500 RCBs exist in our Galaxy. The newly created RGB enriched catalogue is an important step forward 
to significantly increase the number of known RGB stars and therefore better understand their origin. 

Key words. Stars: carbon - AGB and post-AGB - supergiants - circumstellar matter Infrared: stars 



1. Introduction sphere composi tion would help us to constrain simulations of 

merging events (lJeffervetalJl2011l) . 
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C • There exists classes of stars that differ from the vast majority 
^ ' of the stars known in the sense that they are hydrogen-deficient. RGB stars are rare: we currently know about 50 of them 
^ They are called extreme Helium stars (eHe), hydrogen-deficient in the Milky Way (see [Tisserand et al.| |2008l and references 
O Carbon stars (HdC) and R Coronae Borealis (RCB) stars. The therem), which is only a factor of two higher than m the 
^ last two ai-e supergiant carbon-rich stars and have therefore Magellanic Clouds, where 23JlCBs are known dAlcock et al.| 
^ '. strong spectroscopic similarity, but only RCBs are known to un- 2001; Tisserand et alj|2004 12009). The small number of known 
^ • dergo unpredictable fast and large photometric declines (up to 9 ^CB stars and the bias due to their discovery from diff^erent 
> ■ mag over a few weeks) due to carbon clouds formed close to the surveys prevent us from having a clear picture of their true 
• 'H line of sight that obscure the photosphere. Such particular events spatial distribution. Diff^erent views are found in the literature. 
X in such pecuUar stars have made RCBs much-followed objects Iben & Tutukov ( 1985) reported a scale height of h ~ 400 pc 
b ; among many generations of asti'onomers. Nowadays, RCBs have assuming MboI = -5 and concluded that RCBs are part of an old 
■ ri , become even more interesting as they are increasingly suspected disk-like distiibution. However,|Cottrell & Lawson| (|i998|) noted 
to result from the merger of two white dwarfs (one CO- and *at the Hipparcos velocity dispersion of RCB stars is similar to 
one He-), called the double degenerate (DD) scenario. The DD *ose of other cool Hydrogen-deficient carbon stars and extreme 
model has been strongly supported by the obse rvations of an Helium stars, suggesting that RCB stars might have a bulge-like 
'^O overabundance in HdC an d cool RCB stars (IClavton et al.1 distribution. Recendy, adding to the confusion, [Tisserand et al.| 
l2007t iGarcfa-Hernandez et al.ll2010l) . and of surface abundance •I2QQ8) found that the majority of Galactic RCB stars seem to be 
anomalies of a few eleme nts, fluorine in particulai- ( Pandey etiD concentrated in the bulge with the surprising peculiarity of be- 
l2008H.Teff-erv et"d]l20rih . RCBs are actually the favoured candi- '"g distributed in a thin disk structure (61 < h^^^^^^^ < 246 pc, 95% 
dates to be the lower ma ss counterpart of Supernovae type la ob- c.l.). It is therefore necessary to increase the number of known 
jects in a DD scenario dFrver & Diehlll2008t Ibiehl et al. 2008), RCB stars to constrain their spatial distribution and their age, 
therefore detailed studies of their peculiar and disparate atmo- but also understand their past evolution. I note that with an RCB 
phase lifetime of a bout 10 ^ years, as predicted by theoretical 
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per year (iNelemans et al.ll200lt iRuiter et"ani2009h . we can ex- 
pect between 100 and 500 RCB stars to exist in our Gala xy. 

RCB stars are very bright, -5 < My < -3.5 (Tisseran d et al.l 

I2OO9I Fig. 3), and can therefore easily be found anywhere in 
our Galaxy. However, numerous observations are necessary to 
look for their main signature, large declines in luminosity. Well- 
sampled light curves, with a limiting magnitude of ~ 18 for a 
large fraction of the sky, will not be available until the arrival 
of the LSST telescope Q. Fortunately, RGBs are also known to 
possess relatively warm and bright circumstellar shells, easily 
detectable in the mid-infrared. These shells are made of amor- 
phous carbon dust which translates to an almost featureless mid- 
infrared spectrum (Anibal Garcfa-Hernandez et al..2011) , unlike 
the spectra of classical old stars which have siUcate and hy- 
drogen rich dust shells. Therefore, one can imagine finding 
RGBs using only publicl y available infrared ca talogues. Such at- 
tempts have been made : iTisserand et al.l (12011) found two new 
RGBs in the Galactic bulge using mainly mid-in frared SpitzeiE 
GLIMPS^I data and OGLE-Ilfl light curves, and lKraemer et aLl 
(|2005) found cool objects in the Small Magellanic Gloud 
(SMG) with featureless spectra using the Spitzer spectrograph. 
Therefore it is interesting to look at how the first data release 
of the all-sky mid -infrared survey W ISE (Wide-Field Infrared 
Survey Explorer) (I Wright et al.ll20I0l) can improve our search 
for RGBs. This is the goal of the research described in this arti- 
cle. 

In Section |2] I describe briefly the WISE survey and cata- 
logue and explain the impact of the [4.6] band bias observed on 
bright objects using the spectral energy distribution of known 
RGB stars. I describe in Sect.[3]all the criteria used in the analy- 
sis to select a small subsample of the WISE catalogue, enriched 
in RGB stars. Section |4] is a discussion of the outcome of the 
analysis and of the characteristics of the newly formed cata- 
logue. 

2. WISE catalogue and the known RCB stars 

The Wide-field Infrared Survey Explorer (WISE) mapped the 
entire sky during 2010 in 3.4, 4.6, 12 and 22 /im with, re- 
spectively, 6.1, 6.4, 6.5 and 12.0 arcsec in angular resolution 
and 0.08, 0. 11, 1.0 and 6.0 mJy in point source sensitivities 
at 5 sigma (Wright et al.ll2010) . A few months later, a prelim- 
inary data release (WISE-PDRl) was delivered from the first 
105 days of observation. The catalogue contains 257 million 
sources spread over 57% of the sky. Fortunately for an RGB star 
search, the released sky area covers the entire Galactic bulge (see 
Figure[TJ, where most of the known RGBs are located. 

Indeed, 46 of 52 Galactic and 4 of 23 Magellanic RGBs ai-e 
catalogued in WISE-PDRl, but also afl 5 known HdGs and 4 
DY Per type stars. All WISE magnitudes and the 1 -sigma error 
associated are listed in Table |3] for each object. In Table |4] I 
also list fluxes and 1 -sigma errors of all bright Galactic RGBs 
observed by the AKARI satellite, which did a mid-infrared all- 
sky survey in 2006 in 6 bands (centred at 9 and 18 yum with the 
IRG camera, and 65, 9 0, 140 and 160 fim with the FIS camera) 
dMurakami et al.ll2007h . 

' LSST: Large Synoptic Survey Telescope ilvezic et aU2008h 

^ Spitzer is a space telescope launched in 2003 to do infrared imaging 
and spectroscopy (Werner et al. 2004) 

^ GLIMPSE: Galactic Legacy Infrared Mid-Plane Survey 
Extraordinaire (Churchwell et al. 2009) 

'* O GLE-III: Optical Gravitational Lensing Experiment dUdalskil 
120031) 
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Figure 1 : Representation of the Galactic map, b vs /, with all 
objects selected by the analysis (black dots) and the known 
Galactic RGB stars (larger red dots). The black lines represent 
the approximate limits of the sky area released in the WISE- 
PDRl catalogue. 



Photometry of bright non-saturated sources in WISE-PDRl 
has an accuracy of about 2% in [3.4 ], [4.6] and [ 12], and about 
3% for [22] (see for more details Gutriet al.l201lL the WISE pre- 
liminary release explanatory supplement document). The sensi- 
tivity varies significantly due to the different depth of coverage 
(~ 10 epochs on average), the background emission, and the 
source confusion. Saturation begins to affect sources brighter 
than approximately 8.0, 6.7, 3.8 and -0.4 mag respectively in 
all four bands. Most of the known RGBs are brighter than these 
limits, therefore only the PSF-fitting magnitudes will be used in 
the study. Measurements of saturated sources are realised with 
the non-saturated pixels in the objects'wings. Profile-fit pho- 
tometry begins to fail for sources brighter than 1.0, 0.0, -2.0 
and -6.0 mag, which are fortunately brighter than the bright- 
est known RGB star. Photometric bias due to saturation remains 
small (< 0.1 mag) for the [3.4], [12] and [22] WISE bands; how- 
ever, an over-estimate in brightness is observed in the [4.6] band, 
up to nearly 1 mag for objects brighter than 3 mag above the sat- 
uration limit (which corresponds to [4.6] ~ 3.5 mag). The im- 
pact of this bias on bright [4.6] objects will be discussed several 
times in this article, especially in Sect. 12.2.11 

2.1. WISE variability and classification flags 

RGB shell brightness varies. iFeast et al. l (fT997h have shown that 
one can also observe in the L band (3.5 jum) the short-term oscil- 
lations that are observed in the optical (+0.5 mag), as the shell 
reflects light from the photosphere. However, the typical RGB 
large photometric declines are not observed in the mid-infrared 
as only the photosphere gets obscured by clouds. The effect of 
these clouds on the shell luminosity is different. Glouds pro- 
duced in any direction, not just in the line of sight, will obscure 
the light coming from the star and therefore make the shell lu- 
minosity gradually fainter Mid-infrared monitoring could there- 
fore be used as a good indicator of change in the dust formation 
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rate around t he photosphere of an RCB star. A clear example is 
presented bvl Feast et al.l(ll997h with UW Cen, where a variation 
of up to 3 mag in a ~ 100 days was observed. This phenomenon 
is also well described bv lClavton et al.l (Il999h . 

WISE-PDRl gives a variability flag for each band. Therefore 
one can expect to get a positive variability flag for some of 
the RCBs catalogued. However, none of them was flagged to 
be variable, which is surprising as pulsations with periods be- 
tween 30 to 60 days and amplitudes of ~ 0.5 mag should have 
been noticed. The variability flag was calculated if an object 
was observed at least 11 times, but only 45% of the RCBs 
have passed this threshold. (Z Umi has the record with 34 mea- 
surements, followed by R CrB and Y Mus with 20 measure- 
ments.) Furthermore, successive observations of a WISE field 
were made on a short time scale. They span at maximum a time 
range of a few days, which could be too short to observe vari- 
abihty in RCBs. 

A classification (star/galaxy) flag is also available in WISE- 
PDRl, but it will not be used in further analysis as most of the 
catalogued RCBs were found to be not consistent with a point 
source. 

2.2. Spectral energy distributions 

The spectral energy distributions (SEDs) of all bright Galactic 
RCB and HdC stars were reconstructed using catalogues from 6 
different surveys done in the last 10 years (see Figuresl9]tofT2ll. 
In the optical, it is important to measure the maximum bright- 
ness of an RCB and only well sampled monitoring on a long 
time scale can give us a confident result. Therefore, light curves 
ubhshed by AAVSCEI (4 bands : B, V, R and I) and ASAS-3 
(V band) surveys have been analysed. The DENI^ I magni- 
tude was also used if the epoch corresponded to a maximum 
brightness phase, as shown by AAVSO or ASAS-3 lightcurves. 
Overall, the maximum brightness measured in each optical band 
is accurate to ~ 0.05 mag (1-sigma standard error). In the near- 
infrared, the 2MAS^|J, H and K magnitudes were used. Due to 
the pulsating variability of RCB stars brightness, the accuracy of 
the maximum brightness measurement in these 3 bands is low, 
and estimated to be ~ 0.3 mag. A carbon extinction correction 
was applied if the 2MASS epoch corresponded to a fading phase 
of an RCB. This was the case for only 7 stars : SU Tau, UW 
Cen, RZ Nor, V517 Oph, WX CrA, V348 Sgr and DY Per We 
used for such corrections the AV magnitude variation observed 
at that particular epoch and the abso rption coefficients of pure 
amorphous carbon dust presented by iGroenewe gen' ('1995', Fig. 
2). All magnitudes used in the study are listed in Table [T]together 
with the interstellar reddening and extinction factors, E(B - V) 
and A y, obtained from the COBE/DIRB^ map tSchlegel et all 
Il998h with a 4 pixel interpolation. The 1-sigma error on E(B- 
V) is about 0. 1 mag and becomes higher at higher extinction, 
up to 0.3 mag at E(B - V) ~ 1.2 mag (Tisserand etaL.2008.) . 
These extinction coefficients we re app lied to all optical and near- 
infrared magnitudes using IPeil ([T992I) extmction curves, except 



^ American Association of Variable Star Observers, URL: 
|http ://www.aavso.org/vstar/vsots/0100.shtml 

^~A11 Sky Automated Survey jPoimanskil Il997h . URL: 

http://www.astrouw.edu.pl/asas/?page=main 
^ DENIS: The D eep Near-Infrared Southern Sky Survey 

jEpchtein et allll994l) 

* 2MASS: Two Micron All Sky Survey jSkrutskie et alj2006h 

' COBE: COsmic Background Explorer ; DIRBE: Diffuse Infrared 

Background Experiment 



for 7 stars located at low Galactic latitude (|b| < 5 deg) where 
ISchlegel et al.l (Il998h note that the calculated reddening is un- 
certain and untrustworthy as no contaminating sources were re- 
movecE For these stars, the reddening was calculated from the 
Asplund et al.l (1200 0) photosphere effective temperature deter- 
mined using high resolution spectra, and the V - I colour index 
of the stars before reddening correction. 

The SEDs were fitted using the program DUSTY 
(iNenkova et al.l2000l) with two (or three if necessary) black bod- 
ies, a shell m ade entirely of amorphous carbon grains and an 
MRN (Mathis et arill977h size distribution from 0.005 to 0.25 
fim. The DUSTY models used have a photosphere temperature 
ranging from 2800 and 20000 K with a step of 100 K and a 
shell temperature ranging from 300 to 1200 K with a step of 50 
K. The shell visual optical depth, Ty, varied also between 10"^ 
and 10 on a logarithmic scale with 40 increments. The results 
are presented in Figures |9] to [12] together with the fitted effec- 
tive temperatures. (As mentioned earlier, for 7 RCBs located at 
low Galactic latitude, the photosphere e ffective temperature was 
fixed to the temperature determined by Asplun d et al.ir2000h . In 
the mid-infrared, only three WISE magnitudes were used: the 
[4.6] band was not considered. AKARJ magnitudes were fitted 
if no WISE magnitudes were available - this is the case with UV 
Cas, U Aqr and V482 Cyg - or if a clear third cold black body 
was observed (as with DY Cen, MV Sgr, UW Cen and WX CrA). 

As described in Sect 12.11 RCB' shell brightness varies. It 
is therefore not surprising to observe differences in luminosity 
between the WISE and AKARI data, which were taken 4 years 
apart. Such differences are visible on the SEDs presented, partic- 
ularly with RT Nor, where a ~ 2 mag variation is observed, but 
also with V1783 Sgr and V3795 Sgr. I note that the 2MASS K 
band can also be affected by shell brightness variation in the case 
of a hot shell (see WX CrA SED for example). The K flux was 
not used in the SED fit if it clearly contradicted the mid-infrared 
fluxes. 

On average, the accuracy on the temperature and optical 
depth estimates is not better than ~ 10%, because of the com- 
bined effect of the interstellar dust extinction uncertainty (main 
effect), the missing U band magnitude and the fa ct that I did 
not u se a model of a carbon rich photospher^EI (see lLobel et al.l 
11998 !). I compared the photosphere effective temperature found 
to the more reliable one estimated by Asplund et al. (200(1) 
for 17 RCBs using high resolution spectra (temperatures from 
lAsplund et al.l i200 0, are also indicated on the SED Figures) . 
Clearly, the SED method overestimates the photosphere effec- 
tive temperature by about 600 K on average. On the other hand, 
for 4 of the 7 RCB stars (RZ Nor, Y MUs, UV Cas and V3795 
Sgr) located at low Galactic latitude, the best SED model found 
with a photosphere temperature fixed to the one estimated by 
lAsplund et al.l (l2000h and with an interstellar reddening deduced 
from this temperature, does not fit well the near-infrared mag- 
nitudes, indicating that the photosphere temperatures could be 
underestimated. 

Nevertheless, with this simple exercise, I can confirm that 
RCB stars have a wide range of photosphere temperatures, with 
the vast majority having a temperature between ~ 4000 and 8000 
K (spectral type K to F), with a few exceptions being hotter 



Effectively, for 5 of the 7 RCBs, the fitted photosphere effective 
temperatures were higher than 12000 K, when the Schlegel et al. ( 199^ 
map was used to correct for interstellar reddening. This is about 5000 
K higher than expected. 

' ' Carbon molecules CN and C2 create strong absorption features in 
the optical and the near-infrared continuum. 



3 



Tisserand, P.: Tracking down R Coronae Borealis Stars from their mid-Infrared WISE colours 

Table 1: The optical and near-infrared maximum apparent magnitudes determined from the AAVSO, AS AS, DENIS and 2MASS 
datasets. No interstellar extinction was applied to the magnitudes. 
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.26 (0.76)" 


j.o / ^^z.. jj } 


T IW Ten 


10.0 


9.4 


8.8 
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" From Schleael et al. (1998) map, using 4 pixels interpolation. " Calculated using E(B - V)/Av = 3.08 
" Values in parentheses are the reddening E(B - V) and the extinction Ay determined using the photosphere 
Elfective temperatures found by^sglund^^aL (2000). The related stars are close to the Galactic plane (|b| < 5 deg). 



than 10000 K (DY Cen, MV Sgr and V348 Sgr) (iDe Marco etal.1 
I2OOI . RCBs also have a wide range of shell temperatures rang- 
ing from ~ 400 to ~ 1000 K. The typical RCB shell temperature 
is about 700 K with a visual optical depth of Ty ~ 0.4. The dis- 
tributions of these three parameters are presented in Figure |2] 

The SEDs of HdC stars, presented in Figure[T2] confirm that 
their spectral types are F-G, with an effective temperature rang- 
ing from 5500 to 7000 K. As no brightness declines have ever 
been recorded with HdC stars, it is not surprising not to detect a 
mid-infrared excess in their SEDs, but interestingly in the case 
of one HdC, HD 175893, one can fit a second blackbody which 
corresponds to a shell with an effective temperature of ~ 500 
K. HD 175893 could therefore be an RCB star going through a 
long period of low activity, corresponding to a low, if any, dust 
production rate. Furthermore, I note that HD 175893 has passed 
all the RCB selection criteria described below in Sect. |3] 

2.2.1 . WISE [4.6] band bias on bright objects 

An interesting output of this work comes from the difference of 
the [4.6] magnitude published by WISE and the one expected 
from the best SED model found. Figure [3] summarises the situ- 
ation. As mentioned earlier, the WISE team reported a bias for 



bright objects in the [4.6] band, whose brightness could be over- 
estimated by almost 1 mag. I confirm this effect. One can see 
that for RCBs brighter than [4.6] ~ 5.0 mag, there exists a grad- 
ual increase of an overestimation of the [4.6] brightness, up to 
~ 0.9 mag, and a decrease after [4.6] ~ 3.0 mag down to the 
nominal brightness. RCBs whose [4.6] magnitude is certainly 
affected (with 5.0 < [4.6] < 1.5 mag) are indicated in Table [3] 
The WISE team has not yet found the source of this bias. This 
significant effect will be taken into account in the analysis de- 
scribed below. 



3. Analysis 

Using the WISE and 2MASS catalogues, a series of pragmatic 
selection criteria was developed and optimised to create a rea- 
sonable size catalogue enriched in RCB stars located within 50 
kpc. The criteria are described below. The detection efficiency of 
each cut was controlled using 52 known RCBs (see Table|2]i. 

Colour-Colour cuts on all stars detected in every four WISE 
and three 2MASS broadbands are the main selection criteria 
used, as they do not impact on the distance of the RCB. These 
selection criteria are as follows. 
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Figure 2: From left to right, distribution of the fitted effective 
temperature of the RGBs' photosphere and shell, and the shell 
optical depth. 



SED Model [4.6] - WISE [4.6] (mag) 

Figure 3: Diagram representing the measured WISE [4.6] mag- 
nitudes of known RGB stars versus the difference in magnitude 
between the WISE measurement and the best SED model found. 
For objects brighter than [4.6] ~ 5.0, the [4.6] WISE brightness 
is overestimated, reaching ~ 0.9 mag at [4.6] ~ 3.0 mag. The red 
line is a representation of this bias. 



1. The first selection criterion is simply to keep all objects de- 
tected and listed in the WISE catalogue. This process is made 
difficult by the brightness of RGB stars (often saturated) but 
also by the blending effect at low Galactic latitude where 
many RGB stars are expected to be. 52 known RGBs are 
expected to have been catalogued in WISE-PDRl from the 
sky area released (see Figure [1]). However, two are miss- 
ing, OGLE-GG-RGB-1 & -2 which were discovered already 
from their shell brig htness and colours u sing the Spitzer- 
GLIMPSE catalogue dTisserand et alj|201 II) . These 2 RGBs 
are located at very low Galactic latitude (b ~ 3.0 and 1.8 



Figure 4: Diagram of [4.6]-[12] vs [12]-[22]. The black dots rep- 
resent 3% of all objects that have been catalogued in all the 
WISE and 2MASS bands. Larger green dots represent bright 
known Galactic RGB stars (the two outliers are discussed in 
the text); red dots correspond to the Galactic RGBs found 
toward the bulge (with a maiority being located inside the 
bulge) (" Tisserand etal] 120081: IZaniewski e"tal]|2005h : blue dots 
are 4 RGB stars located in the southern-east area of the Large 
Magellanic Gloud. The purple stars correspond to the 5 known 
HdG stars and the circles to DYPer type stars with a colour cod- 
ing identical to RGB stars. The red lines represent the selection 
cuts used in the analysis. 



deg) in a very crowded field. They are clearly bright and dis- 
tinguishable on WISE reference images, in all bands, but the 
WISE source finder did not succeed in retrieving them. This 
first criterion has therefore an important effect on the over- 
all analysis and final outcome completeness. However, this 
issue may be corrected in further WISE data releases. 
The objects selected have to have an entry in all of the four 
WISE bands. The main WISE criterion to accept an object 
in the catalogue is that this object should have been detected 
with a signal-to-noise higher than 7 in at least 1 band and 
detected in 4 frames minimum. At high Galactic latitude, 
the detected source distribution remains relatively uniform 
up to [3.0] < 16.5, [4.6] < 16.0, [12] < 12.5 and [22] < 9.0 
mag. The fact that the faintest Magellanic RGB, EROS2- 
LMG-RGB-5, detected in WISE-PDRl has magnitudes of 
~ 11.9, ~ 11.4, ~ 8.6 and ~ 7.4 mag respectively in all four 
WISE bands, gives us strong confidence that all RGBs lo- 
cated within ~ 50 kpc could be detected. RY Sgr is the only 
RGB star that did not pass this criterion. The problem is re- 
lated to the other end of the brightness source distribution, as 
RY Sgr is the second brightest RGB star known in the opti- 
cal after R GrB, but the brightest in the mid-infrared. RY Sgr 
failed to have an entry in the [12] band. The impact on the 
overall detection efficiency is minimal, as the vast majority 
of the RGBs we are looking for are fainter. 
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3. All selected objects have to have an entry in all J, H and 
K bands of the 2MASS catalogue. All 49 remaining known 
RGBs have passed this selection. The 2MASS point source 
magnitude limits are ~ 15.8, ~ 15.1 and ~ 14.3 mag respec- 
tively in J, H and K. With Mk 6, an RGB star located 

within ~ 50 kpc should have an entry in the 2MASS cata- 
logue if there is not a very strong interstellar extinction on 
the line of sight (A™"^'- < 15 mag). Also, the maximum ex- 
tinction d ue to a cloud o f carbon soot is known to be A ~ 9 
mag in V (lGlavtonll996h . which corresponds to ~ 3 mag in J 
and ~ 1 .6 mag in K. Even in that extreme scenario, an RGB 
star located within ~ 50 kpc should have been detected by 
the 2MASS survey under reasonable interstellar extinction 
conditions (A^™**' <, 4 mag). 

4. The first Golour-Golour selection cut was applied on the 
WISE catalogue to select a reasonable number of objects to 
work with : 



0.75 < [4.6] - [12] < 3.0 & [12] - [22] < 1.3. 



(1) 



The selection is illustrated in Figure |4] One can see that this 
selection is very effective as it selects all objects with a clear 
shell signal and eliminates most of the galaxies located in 
the top-right corner of the diagram. It keeps only 5% of all 
objects catalogued, but rejects only 2 of the RGB stars used 
as references. These 2 RGBs are DY Gen and MV Sgr They 
do not resemble the majority of RGB stars known, as they 
are hot (Tgff ~ 12000 K, see Figure |9]l and surrounded with 
multiple shells. Furthermore, and more interestingly, RGB 
stars are located in an underpopulated area of the diagram. 
They lie in an area well separated from the vast majority of 
the cooler AGB stars. This selection is refined in selection 6, 
and a discussion of the impact of the bias that affects bright 
objects in [4.6] is also presented. 

RGB stars are known to present an excess in t he near-infrared 
compared to classical supergiant stars (see iTisserand et al.l 
12009. Fig.4) due to the warm circumstellar shell that con- 
tributes significantly to the near-infrared fluxes. A selec- 
tion has therefore been applied in the J - H vs. H - K dia- 
gram. Because the interstellar extinction affects these mag- 
nitudes significantly, cuts were developed for the 3 follow- 
ing Galactic latitude ranges. A, B and C: A: b > 2 deg, B: 
1 < b < 2 deg (which corresponds to (Ay) ~ 2 mag) and G: 
b < 1 deg (with (Ay) ~ 5 mag). Three selection cuts were 
applied and are described below. The parameters for each 
Galactic latitude range are given in square brackets as fol- 
lows [A,B,G]: 

(H-K)> [0.2,0.32,0.5] 
if (//-/:)< [0.8,0.92,1.1] : 

{J-H)<{H-K) + [0.2, 0.28, 0.4] 
if {H-K) > [0.8,0.92,1.1] : 

{J-H)<(5mH-K)-\l3 (2) 

The interstellar extinction vector corresponds to 
(J - H)/(H - K) ~ 5/3, therefore the last selection cut 
is identical for each Galactic latitude range. These 3 se- 
lection cuts are represented on Figure |6] Only 2% of the 
remaining objects passed this criterion, and only two RGB 
stars, Y Mus and XX Gam, were eliminated. 
The selection criterion 4 is readjusted to reject most of the 
AGB stars, and the applied cut is illustrated in Figure |6]: 




Figure 5: Golour-magnitude diagram K vs J - K. Black points 
represent the remaining objects after application of selection cri- 
terion 6. The large dots correspond to RGB or HdG stars with the 
same colour coding as in Figure|4] Extinction due to amorphous 
carbon grains (G) and interstellar dust (I) is represented with two 
vectors. The red lines represent selection cut number 7. 



As discussed in Sect. 12.2.11 the brightness of bright objects 
in [4.6] is overestimated. RGBs affected by this bias are in- 
dicated in Table |3] and the estimated corrections are repre- 
sented in Figure |6]with black vertical lines. After correction, 
one can see that none of the RGBs reached the upper limit 
set at 3.0 mag on [4.6] - [12]. The impact of the [4.6] bias on 
this selection criteria is therefore very limited. On the other 
hand, two out of the three RGBs that are eliminated by the 
present selection would have been selected if the bias did not 
exist. As we can expect that many of the RGBs we are look- 
ing for are relatively bright in the [4.6] band, all 3 known 
RGBs rejected will nevertheless be counted in the detection 
efficiency estimation. 

The present criterion target RGB brightness and therefore 
their distance modulus. A selection was applied on the 
colour - magnitude diagram K vs J - K. The colour J - K 
of RGBs is due to three effects: the interstellar reddening, 
the warm circumstellar shell contribution to near-infrared 
fluxes and for some RGBs the reddening due to clouds 
made of carbon soot during an extinction event (at maxi- 
mum A(J - K)cioiid ~ 1.4 mag). The RGBs' J - K colour is 
therefore not simple and the selection limits slope were cho- 
sen with care. They are illustrated on Figure |5] The goal of 
this analysis is to create a catalogue enriched with RGB stars 
located within 50 kpc, therefore all RGBs discovered in the 
Large Magellanic Gloud (LMG) were represented in this dia- 
gram (they all have an entry in the 2MASS catalogue). They 
were used to estimate the lower end of the brightness cut. 
Only 2, EROS2-LMG-RGB-2 & -3, out of 19 LMG RGBs 
were not selected. 



[4.6] - [12] > 1.8([12] - [22]) + 0.03 



(3) 



{5mj-K)+A <K< {5l%){J-K) + \2 



(4) 
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On the bright end of the selection, 2 of the brightest known 
RGBs, R CrB and V854 Cen, are located close to the up- 
per limit. Very bright RCBs are not the primary target of this 
analysis as very few are expected, therefore only a small mar- 
gin was used in this choice. 
8. As mentioned earlier, the RGB's warm circumstellar shell 
can contribute significantly to the near-infrared fluxes. This 
effect can be understood in a J - K vs J - [12] Colour-Colour 
diagram. Effectively, a correlation between both colours 
is observed for RCBs (see Fig. Ell. When the RCB shefl 
becomes brighter, its flux contribution to the K band in- 
creases and as a consequence J - K becomes redder This ef- 
fect is in combination with the interstellar reddening effect. 
Interestingly, it is also relatively simple in that diagram to 
reject the remaining main contaminants, the Mira type stars. 
Indeed, they are located in an isolated area of the diagram. 
The applied cuts were designed to optimise the rejection of 
Miras without impacting too much on the RCB detection ef- 
ficiency. The second selection follows the interstellar red- 
dening vector: 

J -K ^2.25 if (7-[12])<6 (5) 

7-/: < 0.61(7- [12])- 1.43 if {J-[12])>6 (6) 

Only 2 known RCB stars, out of the 42 that passed all the 
precedent criteria, were rejected by the present selection. 
They are EROS2-RCB-CG-12 and MACHO-308-38099.66. 
It is interesting to note that the former one has a light 
curve that presents large oscillations and may therefore be a 
Mira type star, instead of an RCB (see .Tisserand et al...2008k 
Fig.9). 

4. Discussion 

1602 objects have passed the pragmatic selection criteria that 
have just been enumerated. Of the 52 already known RCBs that 
are located inside the sky area covered by the WISE survey pre- 
liminary data release and used as reference, 41 have passed all 
8 criteria, which corresponds to a detection efficiency of about 
77%. For the detection efficiency calculation, it would have been 
preferable to use a Monte-Carlo simulation of RCBs with a range 
of photosphere and shell temperatures, but our knowledge of the 
distribution of these last parameters is very limited. It is worth 
noting that the RCB sample used as reference is biased toward 
bright RCBs and therefore had a negative impact on the detection 
efficiency (see criteria 2 and 6). Bright RCBs are not the primary 
target of this analysis as very few of them, if any, are expected 
to have not been discovered yet. At the other end, the selection 
is biased against RCB stars with faint shells and RCB stars in 
a crowded environment (see criteria 1 and 5). In the last case, 
a detection algorithm correction in a future WISE data release 
may resolve the issue. 

About 70% of the objects selected are located towards the 
Galactic bulge (see Figure[TJ, at +5 deg from the Galactic plane, 
and 85 (~ 5%) are located in the Large Magellanic Cloud, whose 
sky area has been only partially covered in WISE-PDRl. 

The final catalogue was cross-matched with the DENIS and 
USNO-BlFl catalogues, using a strict matching radius of 1 arc- 
sec. Only 895 objects have an entry in the DENIS catalogue and 
538 in the USNO-Bl one. DENIS did a survey of all the south- 
ern sky in I, J and K bands, while USNO-B 1 covers the entire 

'2 USNO: United States Naval Observatory jMonet et al]|2003h 



sky in B, R and I (2 epochs in B and R). Overall, using both 
catalogues, 711 objects have a valid I band magnitude, but sur- 
prisingly only 285 have a valid entry in both catalogues. The I 
band magnitude distribution is presented in Figure [8l as well as 
the difference in magnitude between I peni s and Iusnq-bi- W ith 
an absolute magnitude Mi ~ -5 mag (Tiss erand et al Bool , an 
RCB star located at 50 kpc would have an apparent magnitude 
of 7 ~ 13.5 mag during a maximum brightness phase, which is 
brighter than the median (~ 14) of the distribution presented. The 
Idenis and Iusno-bi magnitudes should be used to prioritise fur- 
ther follow-up. During some large photometric decline phases, 
some RCB stars could have been observed at fainter magnitude 
up to the magnitude limit (iDENis.iim ~ 18.5 mag, or ~ 17 mag 
in crowded area): Figure [8] on the right, shows that some of the 
selected objects have known photometric variation up to 6 mags 
between both survey epochs. 

The fact that only ~ 33% of the 1602 selected objects have 
an entry in the USNO-B 1 catalogue indicates that most of those 
1602 objects are faint, with a visual magnitude lower than 20. 
The interstellar extinction has certainly an inffuence on this low 
percentage, but it is not the absolute answer as half of these faint 
objects are located at 2 degrees or more from the galactic plane 
(\ b \> 2 deg) where extinction Ay is lower than ~ 2 mag. Under 
such extinction, an RCB star at maximum brightness and located 
within 50 kpc is still detectable, therefore these objects could be 
intrinsically optically faint and a contaminant of our research. 
However, they may also be (or a fraction of them) RCB stars that 
are heavily dust enshrouded. If one doesn't expect that during 
the lifetime of an RCB star, there exists a long phase of high 
dust production rate, then these objects should be considered as 
weak candidates. However, it is worth noting that the RCB star 
V854 Cen rem ained 7 mag below its maximum brightness for 
about 50 years (lGrasdalenetal.lll986h . 

I note also that a third of the 1602 selected objects have a 
WISE [4.6] magnitude brighter than 5 mag. The [4.6] brightness 
of these objects is therefore overestimated (see Sect. l2.2.Tl i. 

The selected list of objects was also cross-matched with the 
SIMBAD database with a matching radius of 5 arcsec. Figure|2l 
right side, summarises the situation. 569 of the 1602 objects 
were found to have a classification, but 307 of those are classified 
only as IRAS sources and 30 others only as stars, without further 
informations. However, of the 232 remaining, 62 stars are clas- 
sified interestingly as Carbon stars; they are stars that will need 
to be looked at closely in future photometric and spectroscopic 
follow-ups. Also, 46 objects are classified as Variable, Pulsating, 
Semi-regular or RV Tauri stars. RCB stars are known to present 
periodic pulsating variability with an amplitude of ~ 0.5 mag; 
therefore these stars will be interesting to monitor and study 
more closely. Despite the fact that many of the objects classified 
as Mira or OH/IR stars were rejected using selection criterion 
8, 49 remained in the final sample. Most of them are certainly 
genuine, particularly at J - [12] > 9 mag, but it is also worth in- 
dicating that a few Miras are located in the J - K vs J - [12] di- 
agram at a position not expected for these type of object and can 
therefore be suspected of misclassification. There are 14 such 
objects with J - K < 3 mag. A visual inspection of the ASAS-3 
lightcurves of these objects have shown that 2 of them present 
Mira-type photometric variation (Id 45 and 1 197), 2 remain sta- 
ble (Id 623 and 1409), but more interestingly 3 present clear pho- 
tometric variations typical of RCB stars: Id 683, 793 and 917, 
which are named respectively V653 Sco, 10 Nor and V581 CrA. 
The remaining 7 objects have not been catalogued and followed- 
up by the ASAS-3 survey. I note that 26 objects are classified 
as galaxies or clusters of galaxies. These objects are faint, with 
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Figure 6: In both figures, all points represent objects that have passed the first 4 selection criteria. They are colour coded to represent 
their Galactic latitude b and therefore the interstellar reddening that affects their near-infrared J, H and K magnitudes: black points 
correspond to objects with b > 2 deg, light blue ones are objects with 1 < < 2 deg (corresponding to average extinction of 
(Ay) ~ 2 mag), and purple ones are objects located close to the Galactic plane, with b < I deg (with (Ay) ~ 5 mag). The larger dots 
and stars represent known RGB and HdC stars with identical colour coding to Figure |4] The dashed and solid red lines represent 
the selection cuts number 5 (left) and 6 (right) used in the analysis. Left: J - H vs H - K diagram. Right: [4.6] - [12 vs [12] - [22] 
diagram. The black vertical lines correspond to the brightness correction that needs to be applied to bright RGBs in the [4.6] band 
due to the photometric bias (see Sect. 12.2.11 1. 




Figure?: Both diagrams represent J - K versus J - [12]. Left: the black dots represent the remaining objects after application of 
selection criterion 7. The red line represents the selection cut number 8. The larger dots correspond to RGB or HdG stars that have 
also passed the first 7 criteria, with the same colour coding as in FigureH] The vector represents the interstellar reddening. The black 
curve corresponds to the combination of blackbodies consisting of a 6000 K star and a 700 K shell made of amorphous dust grains, 
in various proportions ranging from all "star" to all "shell". Right: the symbols represent the classification found in SIMBAD. About 
15% of the objects selected with the first seven criteria have an entry in the SIMBAD database with a matching radius of 5 arcsec. 
The dashed line represents the selection cut number 8. 
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Figure 8: Left: distribution of the I magnitude for 711 of 1602 objects selected by the analysis. Only 711 objects have a valid I band 
magnitude in the DENIS and/or USNO-B 1 catalogues. The dotted lines represent the same distribution for known RCB stars. Right: 
Idenis vs Idenis - lusNO-Bi for the 285 of 1602 objects that have a valid I magnitude in both catalogues. The difference in magnitude 
between both catalogues covers about 9 magnitudes, showing that many selected objects undergo large photometric variation. 

Table 2: Number of selected objects after each selection criterion. 



Selection criterion 


Number of 


Number of 


RCB stars 




WISE objects 


known RCBs 


eliminated 


0: Located in WISE-PDRI sky area 




52 




I: Catalogued by WISE 


257310278 


50 


OGLE-GC-RCB-1 &-2 


Catalogued in the [22] band 


12192351 


50 




2: Catalogued in all four WISE bands 


9636678 


49 


RYSgr 


3: Catalogued in all three 2MASS bands 


6832051 


49 




4: Cut on ([4.6] - [12] vs [12] - [22]) 


337359 


47 


DY Cen and MV Sgr 


5: Cut on (I - H vs H - K) 


9660 


45 


Y Mus and XX Cam 


6: Cut on ([4.6] - [12] vs [12] - [22]) 


4146 


42 


EROS2-RCB-CG-3 &l -5 and SV Sge 


7: Cut on (K vs J - K) 


3058 


42 




8: Cut on (J - K vs J - [12]) 


1643 


40 


EROS2-RCB-CG-12 and MACHO-308-38099.66 


Final (without the known RCBs + HD 175893) 


1602 







K > \2 mag, indicating that galaxies could therefore be a major 
contaminant below this limit (270 on the 1602 selected objects 
are fainter than K > 12). Finally, 17 objects are classified as 
Emission stars; such objects could be contaminants in our search 
for RCB stars, but it is interesting to note that hot RCBs, such 
as DY Cen, MV Sgr and V 348 Sgr, present also an emission 
type spectrum (iDe Marco et al. 2002 ). Furthermore, during faint 
phases, RCBs spectrum generally presents many emission lines 
(IClavtonlll996h . 

There does not exist any reference sample of stars to esti- 
mate the fraction of RCB stars that are potentially in the cata- 
logue. Furthermore, little is known on the different phases during 
the lifetime of an RCB star and as mentioned earlier, a signifi- 
cant fraction of RCB stars could be heavily dust-enshrouded and 
therefore faint. However one can use the objects that have been 
classified in the SIMBAD database to obtain a crude estimate of 
this fraction: ~ 15% of these objects are RCBs, it would then 
corresponds to about 240 new RCBs. 

The RCB enriched catalogue is available from the following 
URL: http://www.mso.anu.edu.au/~tisseran/RCB/ and will also 
be available through the VizieR0 catalogue service. A short ver- 
sion of it is given as example in Table |5] The equatorial and 
Galactic coordinates as well as all four WISE, three 2MASS, 



URL: http://vizier.u-strasbg.fr/viz-bin/VizieR 



three DENIS magnitudes and their 1-sigma errors are listed for 
all 1602 objects selected. The five USNO-B 1 magnitudes are 
also listed, but not the individual measurem ent error as they w ere 
not delivered in the original catalogue (see iMonet et al.l2003L for 
an estimate of the photometric accuracy). If one magnitude was 
not available, its value was replaced by the number -99. Also, if 
more than one epoch were available in the DENIS or USNO-B 1 
catalogues for one particular object, only the epoch related to the 
brightest magnitudes was kept. The last column of the catalogue 
gives the SIMBAD classificatioiQ as of July 201 1, found using 
a 5 arcsec matching radius. An underscore character was given 
to the objects that had no classification in SIMBAD. 

5. Conclusion 

Using both the 2MASS and WISE preliminary data release cat- 
alogues, I have created a catalogue enriched with Galactic RCB 
stars that lie within a distance ~ 50 kpc. The selection criteria to 
select a subsample of objects with mid- and near- infrared prop- 
erties similar to RCBs, were mostly based on colour-colour di- 
agrams. This catalogue contains only 1602 entries, with ~ 70% 
of the selected objects being located toward the Galactic bulge. 



URL of the object types in SIMBAD (classification version : 12- 
Jul-201 1): http://simbad.u-strasbg.fr/simbad/sim-display ?data=otypes 
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A reference sample of 52 known RCB stars was used to moni- 
tor the detection efficiency of the selection. About 77% of them 
were recovered. Such a high detection efficiency gives strong 
support to the RCB content of this catalogue. Each selected ob- 
ject will now need to be followed-up spectroscopically to dis- 
cover its true nature. It is encouraging to observe that 3 of these 
1602 selected objects. Id 683, 793 and 917 (resp. V653 Sco, lO 
Nor and V581 CrA) were found to present ASAS-3 lightcurves 
with brightness variations typical of RCB stars , but are misclas- 
sified i n the SIMBAD database as Miras. Also, Kiibunchoo et al.l 
(1201 lb has recently discovered that NSV 11154 is a new RCB 
stars. This star is listed in the catalogue under the Id 1240. 

An analysis of the Spectral Energy Distribution of known 
RCB stars confirms that RCBs effective temperatures range 
mostly between 4000 and 8000 K (spectral type K to F), with 
a few exceptions being hotter than 10000 K. The RCB shell ef- 
fective temperatures range between 400 and 1000 K. The typical 
shell temperature being ~ 700 K with a visual optical depth of 
TV ~ 0.4. HD 175893 is the only HdC star (spectral type G to 
F) that presents an infrared excess indicating the existence of a 
warm circumstellar shell. Furthermore, HD 175893 has passed 
all the selection criteria to create the RCB enriched catalogue, 
therefore HD 175893 should be considered to be an RCB star 
under a phase of low activity in terms of dust production. 
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Table 3: WISE magnitudes of catalogued RCB, HdC and DY Per stars. 



Name 


[3.0] 


O" [3.0] 


[4.6] 


""[4.6] 


[12] 


C^[12] 


[22] 


0"[22] 


Galactic RCB stars 


ui Cen 


1 n y1 Q o 
I0.4J8 


0.0/3 


9.1 /6 


O.Ozl 


4.1U6 


0.014 


2.39s 


0.013 


ES Aql 


5.954 


0.044 


4.660 


0.035 


3.344 


0.020 


2.837 


0.017 


FH Set 


6.076 


0.044 


4.81 1* 


0.031 


3.589 


0.019 


3.004 


0.021 


GU Sgr 


5.993 


0.051 


4.289* 


0.039 


3.021 


0.019 


2.624 


0.016 


MV Sgr 


8.144 


0.024 


7.400 


0.021 


4.695 


0.019 


2.234 


0.017 


R CrB 


3.505 


0.067 


1.920* 


0.002 


-0.273 


0.003 


-0.833 


0.019 


RS Tel 


7.075 


0.031 


5.672 


0.027 


3.676 


0.018 


3.057 


0.021 


RT Nor 


5.877 


0.044 


3.987° 


0.035 


2.856 


0.020 


2.628 


0.018 


RZ Nor 


6.687 


0.030 


4.820* 


0.031 


2.836 


0.021 


2.087 


0.018 


RY Sgr 


2.615 


0.004 


1.434* 


0.001 


— 


— 


-0.564 


0.012 


S Aps 


6.064 


0.044 


5.498 


0.027 


3.757 


0.018 


2.854 


0.018 


SU Tau 


4.928 


0.055 


3.041* 


0.028 


1.484 


0.008 


0.987 


0.013 


SV Sge 


5.600 


0.048 


5.025 


0.033 


3.322 


0.019 


2.191 


0.015 


UW Cen 


5.705 


0.029 


3.322* 


0.016 


1.442 


0.006 


0.655 


0.012 


V1157 Sgr 


6.291 


0.038 


4.582* 


0.031 


3.118 


0.022 


2.549 


0.020 


V1783 Sgr 


6.128 


0.041 


5.082 


0.027 


3.133 


0.021 


2.333 


0.018 


V2552 Opn 


7.360 


0.027 


6.224 


0.023 


4.760 


0.019 


4.285 


0.028 


V348 Sgr 


5.406 


0.020 


3.042* 


0.009 


1.838 


0.013 


1.361 


0.015 


V3795 Sgr 


7.153 


0.029 


5.439 


0.025 


3.045 


0.020 


2.410 


0.017 


V4017 Sgr 


7.276 


0.028 


6.345 


0.023 


4.653 


0.018 


3.941 


0.021 


V517 Oph 


4.969 


0.058 


3.121* 


0.033 


2.031 


0.018 


1.624 


0.016 


V739 Sgr 


6.448 


0.037 


5.283 


0.028 


3.818 


0.019 


3.225 


0.021 


V854 Cen 


3.120 


0.011 


1.813* 


0.002 


0.448 


0.003 


0.398 


0.011 


V CrA 


6.008 


0.050 


4.211* 


0.046 


2.387 


0.012 


1.624 


0.013 


VZ Sgr 


7.405 


0.028 


6.474 


0.019 


4.938 


0.019 


4.260 


0.027 


WX CrA 


6.206 


0.041 


5.100 


0.031 


3.798 


0.017 


3.259 


0.023 


XX Cam 


5.431 


0.053 


4.760* 


0.033 


3.427 


0.019 


2.867 


0.018 


YMus 


8.133 


0.025 


7.855 


0.022 


5.846 


0.018 


4.612 


0.025 


ZUmi 


5.982 


0.045 


4.673* 


0.030 


3.429 


0.018 


2.933 


0.017 


MACHO- 135. 27 132.51 


8.337 


0.024 


7.033 


0.019 


5.043 


0.013 


4.323 


0.031 


MACHO-301.45783.9 


8.540 


0.027 


7.177 


0.021 


5.351 


0.018 


4.578 


0.030 


MACHO-308. 38099. 66 


7.647 


0.023 


6.495 


0.021 


4.829 


0.016 


4.071 


0.025 


MACHO-401.48170.2237 


5.901 


0.044 


4.576* 


0.035 


3.407 


0.019 


2.858 


0.028 


EROS2-CG-RCB-1 


6.666 


0.039 


5.036 


0.026 


3.216 


0.019 


2.334 


0.018 


EROS2-CG-RCB-3 


5.768 


0.047 


3.954* 


0.068 


3.099 


0.018 


2.585 


0.020 


EROS2-CG-RCB-4 


6.41 1 


0.036 


5.080 


0.030 


3.635 


0.017 


2.998 


0.027 


EROS2-CG-RCB-5 


6.319 


0.035 


4.867* 


0.030 


3.784 


0.015 


3.144 


0.022 


EROS2-CG-RCB-6 


7.140 


0.031 


5.933 


0.025 


4.314 


0.020 


3.670 


0.034 


EROS2-CG-RCB-7 


7.210 


0.031 


6.143 


0.023 


4.628 


0.018 


3.985 


0.029 


EROS2-CG-RCB-8 


7.050 


0.032 


5.795 


0.026 


4.189 


0.019 


3.595 


0.036 


EROS2-CG-RCB-9 


7.250 


0.028 


5.406 


0.027 


3.423 


0.018 


2.649 


0.023 


EROS2-CG-RCB-10 


6.739 


0.036 


4.485* 


0.039 


2.546 


0.023 


1.958 


0.022 


EROS2-CG-RCB-11 


6.581 


0.035 


5.456 


0.026 


4.138 


0.016 


3.523 


0.024 


EROS2-CG-RCB-12 


8.067 


0.026 


7.380 


0.022 


6.452 


0.019 


6.263 


0.055 


EROS2-CG-RCB-13 


6.839 


0.034 


5.541 


0.026 


3.985 


0.016 


3.378 


0.026 


EROS2-CG-RCB-14 


6.375 


0.034 


4.695* 


0.031 


3.292 


0.019 


2.678 


0.019 


Magellanic RCB stars 


MACHO- 12. 10803.56 


10.287 


0.025 


9.367 


0.021 


7.744 


0.020 


7.208 


0.055 


EROS2-LMC-RCB-4 


10.001 


0.025 


8.551 


0.021 


6.637 


0.018 


6.139 


0.034 


EROS2-LMC-RCB-5 


11.854 


0.025 


11.365 


0.022 


8.623 


0.022 


7.437 


0.067 


EROS2-LMC-RCB-6 


10.604 


0.025 


9.082 


0.022 


6.857 


0.018 


6.102 


0.035 


HdC stars 


HD137613 


5.090 


0.059 


4.909* 


0.033 


5.000 


0.018 


4.854 


0.028 


HD 148839 


6.787 


0.033 


6.804 


0.023 


6.731 


0.019 


6.770 


0.063 


HD 173409 


8.095 


0.018 


8.067 


0.019 


8.021 


0.025 


8.508 


0.349 


HD 175 893 


7.165 


0.029 


6.832 


0.021 


5.104 


0.019 


4.175 


0.025 


HD 182040 


4.953 


0.054 


4.511* 


0.033 


4.730 


0.019 


4.739 


0.032 


DY Per stars 


DYPer 


3.030 


0.089 


2.254* 


0.007 


1.941 


0.009 


2.221 


0.022 


EROS2-CG-RCB-2 


8.189 


0.035 


7.856 


0.030 


7.388 


0.045 


7.376 


null 


MACHO- 15. 10675. 10 


9.525 


0.023 


9.344 


0.020 


8.751 


0.021 


8.679 


0.135 


EROS2-LMC-DYPer-5 


10.304 


0.024 


10.178 


0.021 


9.616 


0.027 


9.091 


0.240 



*[4.6] magnitude affected by bias (see Fig.O 
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Figure 9: Spectral energy distributions of known bright Galactic RGB stars, normalised to flux in V. Red dots represent fluxes in the 
optical (B, V, R and I) and the near-infrared (J, H and K) obtained from AAVSO and the ASAS, DENIS and 2MASS surveys. The 
blue and black dots represent mid-infrared fluxes from the WISE and AKARl surveys respectively. The red line is the best fit found 
with DUSTY models (see text for more details); the related eff'ective temperatures found (~ 10% level accuracy) are listed at the 
bottom left corner in the following order, from top to bottom : photosphere, first and second circumstellar shell. If the photosphere 
effective temperature is bet ween square brack ets, its value was fixed during the fit and corresponds to the photosphere eff'ective 
temperature determined by Asplund et al.l (12000 ). indicated below the name for 17 RGB stars. The values in brackets on the right 
side of the shell temperatures correspond to the visual optical depth found. The broken black line represents a simple blackbody 
function with the photosphere temperature. 
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Figure 10: Spectral energy distributions of known bright Galactic RGB stars, normalised to flux in V. Same caption as Figure|9] 
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Figure 11: Spectral energy distributions of known bright Galactic RGB stars, normalised to flux in V. Same caption as Figure|9] 



14 



Tisserand, P.: Tracking down R Coronae Borealis Stars from their mid-Infrared WISE colours 
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Figure 12: Spectral energy distributions of the 5 HdC stars known, normalised to flux in V. Red dots represent fluxes in the optical 
and the near-infrared obtained from AAVSO and the ASAS, DENIS and 2MASS surveys. The blue and black dots represent the 
mid-infrared fluxes from the WISE and AKARI surveys respectively. The red line is simply a spline function that connects the 
different fluxes. The broken black line represents a simple blackbody function with the temperature of the photosphere. In the case 
of HD175893, the dotted line represents the sum of two blackbodies, the photosphere and the shell. The black bodies' eff'ective 
temperatures of the photosphere (top) and the circumstellar shell (bottom) are listed at the bottom left corner. 

Table 4: AKARI fluxes (mJy) of catalogued bright Galactic RGB stars and DY Persei. 



Name 


[9.0] 


C [9.0] 


[18] 


""1181 


[65] 




[90] 


C [90] 


[140] 


O" [140] 


[160] 


C [160] 


DYCen 


0.538 


0.022 


0.889 


0.016 


















ES Aql 


1.935 


0.127 


0.875 


0.036 


















FHSct 


1.397 


0.032 


0.661 


0.053 


















GUSgr 


1.290 


0.140 


0.751 


0.100 


















MV Sgr 


0.330 


0.019 


1.008 


0.008 






0.496 


0.103 










RCrB 


52.990 


2.440 


21.480 


0.029 


2.656 


0.170 


1.494 


0.114 






0.606 


0.333 


RSTel 


1.870 


0.102 


0.767 


0.011 


















RTNor 


0.257 


0.037 


0.228 


0.021 


















RZNor 


3.119 


0.252 


1.695 


0.076 


















RYSgr 


48.000 


3.660 


20.180 


1.020 


3.495 


0.104 


2.605 


0.139 










S Aps 


1.971 


0.016 


1.001 


0.021 


















SUTau 


14.720 


0.050 


6.161 


0.046 


0.351 




1.179 


0.080 










SV Sge 


3.849 


0.441 


2.293 


0.095 


















U Aqr 


1.260 


0.142 


0.750 


0.002 


















UV Cas 


1.125 


0.046 


0.802 


0.026 


















UWCen 


9.760 


0.070 


5.697 


0.047 


6.645 


0.413 


7.302 


0.332 


4.179 


0.349 


2.809 


1.020 


UX Ant 


0.101 


0.007 






















V1157 Sgr 


2.569 


0.004 


1.145 


0.015 


















V1783 Sgr 


3.918 


0.252 


1.900 


0.032 


















V348 Sgr 






2.813 


0.028 






1.638 


0.018 










V3795 Sgr 


4.574 


0.144 


1.951 


0.069 


















V482 Cyg 


1.128 


0.037 


0.459 


0.020 


















V517 Oph 


6.756 


0.243 






















V739 Sgr 


1.638 


0.065 


0.840 


0.001 


















V854 Cen 


22.970 


1.170 


7.364 


0.033 






0.705 


0.036 










VCrA 


3.606 


0.210 


2.166 


0.013 


















VZ Sgr 






0.283 


0.098 


















WXCrA 


1.742 


0.097 


0.810 


0.001 






1.548 


0.105 


2.919 


6.340 






XX Cam 






1.391 


0.065 


















YMus 


0.230 


0.017 






















ZUmi 


1.724 


0.029 


0.765 


0.024 


















DYPer 


12.210 


0.823 


2.600 


0.157 
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Table 5: First 10 rows of the published catalogue. 



ID 


Ra 


Dec 


1 


b 


WISE 


WISE 


WISE 


WISE 


2MASS 




(deg) 


(deg) 


(deg) 


(deg) 


[3.0] 


f [3.0] 


[4.6] 


f [4.6] 


[12] 


0"[12] 


[22] 


O" [22] 


J mag 


o-j 


1 


105.4079971 


-41.4713669 


251.85892 


-15.80245 


4.852 


0.058 


3.017 


0.016 


1.736 


0.017 


1.197 


0.016 


11.324 


0.024 


2 


259.1376648 


-33.8375053 


352.49854 


2.49227 


8.529 


0.029 


6.471 


0.022 


3.917 


0.016 


2.809 


0.019 


17.348 


-99 


3 


109.6635208 


-24.9482384 


238.16979 


-5.55101 


9.701 


0.026 


8.951 


0.021 


7.202 


0.020 


6.375 


0.085 


11.004 


0.022 


4 


107.7039337 


-6.3668070 


220.80645 


1.38275 


6.906 


0.029 


4.715 


0.034 


2.462 


0.022 


1.557 


0.021 


13.736 


0.030 


5 


273.2704468 


-25.0299644 


6.43576 


-3.35900 


10.990 


0.078 


11.098 


0.085 


9.010 


0.158 


8.046 


0.329 


12.274 


-99 


6 


258.5602112 


-21.4371471 


2.42136 


10.03461 


6.800 


0.032 


5.829 


0.022 


4.418 


0.019 


4.017 


0.023 


9.783 


0.024 


7 


273.8857422 


-25.5371208 


6.25302 


-4.08940 


10.657 


0.050 


10.018 


0.040 


8.367 


0.031 


7.509 


0.178 


12.799 


-99 


8 


273.0597229 


-25.4018135 


6.01729 


-3.36891 


10.507 


0.068 


9.501 


0.038 


7.216 


0.025 


6.009 


0.070 


12.616 


0.029 


9 


272.6348572 


-25.3729763 


5.85848 


-3.01816 


9.329 


0.035 


9.041 


0.026 


8.114 


0.030 


8.043 


0.289 


10.410 


0.021 


10 


273.3207092 


-25.5264320 


6.02011 


-3.63545 


9.215 


0.028 


7.982 


0.021 


6.220 


0.018 


5.708 


0.043 


14.477 


0.068 



Colunms continued 





2MASS 


2MASS 




DENIS 




DENIS 


DENIS 






USNO-Bl 






SIMBAD 


H mag 




Kmag 


(Tk 


I mag 




J mag 


o-j 


Kmag 


Tk 


Bl mag 


B2 mag 


Rl mag 


R2 mag 


I mag 


Class. 


9.154 


0.024 


7.287 


0.027 


15.768 


0.06 


12.201 


0.06 


8.076 


0.06 


-99 


17.06 


-99 


16.97 


15.91 


c* 


16.261 


-99 


13.424 


0.050 


-99 


-99 


-99 


-99 


12.049 


0.11 


-99 


-99 


-99 


-99 


-99 




10.798 


0.026 


10.549 


0.023 


-99 


-99 


-99 


-99 


-99 


-99 


-99 


-99 


-99 


-99 


-99 


*iC 


10.807 


0.027 


8.466 


0.023 


-99 


-99 


13.650 


0.09 


8.349 


0.07 


-99 


-99 


-99 


-99 


-99 




11.694 


0.049 


11.312 


0.038 


-99 


-99 


-99 


-99 


11.294 


0.08 


-99 


-99 


-99 


-99 


-99 




9.036 


0.024 


8.479 


0.021 


10.806 


0.04 


9.393 


0.07 


8.261 


0.25 


14.13 


11.90 


15.26 


11.89 


10.42 


c* 


12.844 


0.064 


12.196 


0.034 


-99 


-99 


-99 


-99 


-99 


-99 


-99 


-99 


-99 


-99 


-99 




12.018 


0.027 


11.493 


0.021 


14.293 


0.04 


12.961 


0.06 


11.776 


0.09 


16.67 


14.12 


16.13 


15.39 


13.81 




10.008 


0.025 


9.735 


0.021 


12.816 


0.32 


11.404 


0.40 


10.819 


0.55 


14.20 


12.67 


14.15 


13.05 


13.38 




12.586 


0.094 


11.067 


0.039 


16.327 


0.09 


13.232 


0.13 


10.515 


0.13 


-99 


-99 


-99 


-99 


-99 


V* 



